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ABSTRACT 

The Schrodinger-Newton [SN] equation describes the effect of self-gravity on the evolution of a 
quantum system, and it has been proposed that gravitationally induced decoherence drives the 
system to one of the stationary solutions of the SN equation. However, the equation by itself 
lacks a decoherence mechanism, because it does not possess any stochastic feature. In the present 
work we derive a stochastic modihcation of the Schrodinger-Newton equation, starting from the 
Einstein-Langevin equation in the theory of stochastic semiclassical gravity. We specialize this 
equation to the case of a single massive point particle, and by using Karolyhazy’s phase variance 
method, we derive the Diosi - Penrose criterion for the decoherence time. We obtain a (nonlinear) 
master equation corresponding to this stochastic SN equation. This equation is however linear at 
the level of the approximation we use to prove decoherence; hence the no-signalling requirement 
is met. Lastly, we use physical arguments to obtain expressions for the decoherence length of 
extended objects. 


I. INTRODUCTION 


A problem of long-standing interest has been to understand if gravity can help solve the quantum 
measurement problem; or even if it does not actually cause collapse of the wave-function, can 
it at least be a source of decoherence. The physical picture being that the gravitational field 
of space-time cannot be exactly classical but must possess intrinsic quantum uncertainty. Can 
this uncertainty relate to the quantum evolution described by the Schrodinger equation in such 
a way that for macroscopic masses the space-time uncertainty induces decoherence and position 
localisation? The answer seems to be in the affirmative, as has been demonstrated in various studies 
by Karolyhazy [l| and then followed by Karolyhazy and collaborators in several papers P-[ll|. and 
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also shown by Diosi and collaborators 


12l-ll5l|. These authors demonstrate, subject to crucial 


assumptions about how quantum matter affects gravity, that there is an inherent uncertainty in 
the space-time geometry, whose origin lies in the quantum nature of the sources which produce 
them. This uncertainty is modelled by introducing a classical stochastic potential in the Schrodinger 
equation, and it is shown that the stochasticity induces decoherence, whose properties depend on 
the mass and size of the quantum object under study. The work of Karolyhazy and of Diosi has 
recently been compared in [l^ . 

Having considered the effect of the space-time uncertainty produced by other objects, on the 
Schrodinger equation, it is also important to consider the effect of self-gravity of a quantum object 
on its Schrodinger evolution. This self-gravity has an intrinsic quantum uncertainty and although 
we do not quite know at present how to describe the self-gravity of a quantum object, one can 
attempt to model it. One well-known approach is the Schrodinger-Newton [SN] equation. It is 
proposed that the self-gravitational potential V produced by a quantum source in state 'k satisfies 
a semiclassical Poisson equation 

VV = 47rGm|Tp (1) 


whose solution is incorporated in the potential dependent Schrodinger equation 

ih—— =--|- mHT 

dt 2m 


This gives the Schrodinger-Newton equation 


13-l2C| 






|T(r',t)|^ 


dt 'Zm J \r — r' 

This equation has been discussed extensively in the literature 


■(ir'T(r,t) 


( 2 ) 


( 3 ) 


21I431I|. The equation is nonlinear 


and deterministic. This makes it different in character from the Karolyhazy equation and the Diosi 
equation, where by virtue of the stochastic potential the models predict decoherence at the level 
of the master equation. The SN equation does not predict decoherence but is instead suggested 
as an equation whose stationary solutions are the ones to which the system will proceed upon 
decoherence, once a gravity-based decoherence mechanism has been incorporated in this system. 
What has been shown is that there is a gravitationally induced inhibition of dispersion of an 


expanding wave-packet, at the critical length Oc ~ h‘^/Gm^ 


22l |. Another concern which can be 


raised about the SN equation is that its deterministic non-linear structure allows for superluminal 
signalling. For a recent critique of the SN equation see 291 and Q- 

The reason why decoherence is not observed in the SN equation is evident: one is only taking 
into account the mean potential in the semiclassical Poisson equation, whereas the Karolyhazy 
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and Diosi models incorporate stochastic fluctuations by way of the stochastic potential in the 
Schrodinger equation, which is modelled after the space-time geodesic uncertainty. We need to take 
into consideration the stochastic fluctuations, apart from the mean, and modify the SN equation 
accordingly. One possible way to do this is to start from the theory of stochastic gravity 
takes into account corrections to the semiclassical Einstein equations 


32|. which 


1 SttG 


(4) 


by considering the role of the two-point fluctuations of the energy-momentum tensor [going from 
semiclassical gravity to the Einstein-Langevin equation]. The purpose of the present paper is to 
develop this idea and propose a new stochastic SN equation. The effect of quantum fluctuations 
of the energy-momentum tensor on the geometry is modelled by defining a classical stochastic 
field. This stochastic field eventually acts as a source for the gravitational potential (which is now 
stochastic) via a modified semiclassical Poisson equation. This stochastic potential is included in 
the Schrodinger equation, in the spirit of Karolyhazy’s and Diosi’s work, and as we shall show, it 
leads to a stochastic SN equation which produces decoherence in the quantum evolution. 

The plan of the paper is as follows. In Section II we recall the Einstein-Langevin equation. In 
Section III we use the Einstein-Langevin equation to derive a new stochastic SN equation. Section 
IV uses this new equation to provide a proof of gravitational decoherence for the gaussian state of 
a point particle. In Section V we obtain elementary estimates for decoherence scales for extended 
objects while working with the SN equation. Concluding remarks are in Section VI. Appendix I 
derives the master equation for the density matrix, corresponding to the stochastic SN equation. 
Appendix II gives details of the decoherence calculation reported in Section IV. 

A stochastic modification of the Schrodinger equation has also been proposed in 
contrast our work with that of 33|| in the concluding section. The application of stochastic gravity 


3^. We 


to the measurement problem has recently been studied in 


34l |. The possible role of gravity in the 


decoherence / collapse of the wave function has been reviewed, amongst other places, in 


35|-l37|. 


II. SEMICLASSICAL GRAVITY AND THE EINSTEIN-LANGEVIN EQUATION 

Semiclassical gravity describes the interaction of the gravitational field, which is treated classi¬ 
cally, with quantum matter fields. The field equation for the classical metric is the semiclassical 
Einstein equation, which gives the back reaction of the matter fields on the spacetime; it is a 
generalization of the Einstein equation where the source is the expectation value in some quantum 
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state of the matter stress-energy tensor operator. The coupling of gravitational field to matter is 
modelled by the semiclassical Einstein equations (j^. 

We will be interested in the weak-field limit 


5 /^ 1 / — T hfj,i/, 


(5) 


where \hf^^\ <C 1 and is the Minkowski metric; thus we obtain the linearised theory of gravity. 
In a weak field situation one can expand the field equations in powers of h^^ using a coordinate 
frame where Eqn. ([5|) holds and one can retain only linear terms. Using Eqns. (j3|) and ([5]), one 


can arrive at the differential equation for the metric field given by 

IQttG 


301] 


G\h^i/ — 




( 6 ) 


In the Newtonian limit the ('kjrool'I') component is large compared to other components of the 
stress-energy tensor. Hence in this limit the field equation for linearised gravity reduces to a 
semiclassical Poisson equation of the form 

AttG 


VW = 


-('kjTool'h), 


(7) 


where /ioom/2 = —V. Using Too = cfp we get the familiar potential held used in the Schrodinger- 


Newton equation, as explained in 


30l |. The same result also follows by noting that for a single 


particle we have g = m|r)(r|. In this manner we can think of the SN equation ([3]) as the non- 
relativistic limit of semiclassical gravity. 

The semiclassical Einstein equations have as a source only the ‘mean held’ coming from the 
expectation value of the quantum stress-energy tensor. It has been argued that there can be 
circumstances where quantum huctuations of the stress-energy tensor about the mean can be 


important, and then the semiclassical theory should be modihed. One possib 


so-called theory of stochastic gravity. Here, the Einstein-Langevin equation [321. l3S 


e modihcation is the 


39| | includes also 


the back reaction on the space-time metric of the lowest order stress-energy quantum huctuations. 
This results in an effective theory which predicts linear stochastic corrections to the semiclassical 
metric. The equations include a Gaussian stochastic tensor held representing the stress- 

energy huctuations, and having the following three properties: 

1. The stochastic average (.}s of the held vanishes i.e. 


— 0 . 


( 8 ) 
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2. The two point correlation of this field is given by 

s ^afifii/{x, y} ■ (9) 

Here x and y are space-time four vectors. Nayij,u{x,y) is called the noise kernel. It is related 
to the stress-energy tensor in the following way 

= {{iay{x),i^u{y)}), (10) 

where {.} is the anti-commutator, (.) is the expectation value and tap = Tap — {Tap)- 

3. Higher moments and cumulants of this stochastic field vanish. 

The metric continues to be classical, and now becomes a stochastic tensor field, sourced by the 
mean as well as by the fluctuations. The Einstein-Langevin equation takes the form 

+ ( 11 ) 

with the last term on the right hand side incorporating the effect of matter fluctuations. 

We will now see how the inclusion of matter fluctuations modifies the SN equation, leading to 
our proposal for a new stochastic SN equation. 


III. A PROPOSAL FOR A STOCHASTIC SCHRODINGER-NEWTON EQUATION 


We consider the weak field limit of the Einstein-Langevin equation (HH). On applying the 
weak-field condition Eqn. ([5|) in (Hip , we obtain the following field equation for the metric 




In the Newtonian limit, we get the equation 

AttG 


V^V{r,t) = 


1 


(Tlfbol'k) -h 2(^00 + aiapy^"^) 


( 12 ) 


(13) 


c" \.- 2 

Now if we consider only the ^oo component (since Newtonian limit has been assumed), we get 

VV(r,t) = ^((T|roo|'k)+eoo). (14) 

As expected, we now have a stochastic source contributing towards the gravitational potential. 
Using Green’s function, we derive the modified potential which includes a stochastic component 


= -Gm f MWT*' - L [ I 22 T 4 )*'. 


r — r' 


r — r' 


(15) 
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Using this potential in the Schrodinger equation ([2]) gives the new stochastic Schrodinger-Newton 
equation 






dr''I'(r, t) — 


Gm f ^oo{r',t) 


/ 


r — r 


t) (16) 


|r —r'l 

This equation can also be worked out by going from field theory to quantum mechanics of N 
particles in non-relativistic limit and then obtaining the SN equation for single particle, following 
the method outlined in 30(]. The modified potential, where the modifying stochastic component 


originates in the fluctuations of the stress tensor, is responsible for the stochastic SN equation. 

A couple of striking similarities with Diosi’s model and Karolyhazy’s model can be immediately 
recognised. First, the total stochastic potential V{r,t) is a linear sum of the Newtonian potential 
and stochastic noise (which is Gaussian). Second, the average of the stochastic component of the 
potential vanishes due to the first property of ^^v{x). In a sense then, the stochastic SN equation 
(unlike the SN equation itself) is on the same footing as the stochastic Schrodinger equations of the 
Karolyhazy model and the Diosi model. The key difference between the first and the latter two is 
that while one accounts for effects of self-gravity, the other two account for stochastic fluctuations 
induced by extrinsic space-time uncertainty. To our understanding, a more general treatment 
should take into account effects of both self-gravity and extrinsic uncertainty in the same scheme. 
We hope to pursue such an investigation in future work. 

For use in the next section, we work out the two point correlation of the stochastic potential. 
Let the stochastic component of the gravitational potential in Eqn. ()15p be denoted byI4i(r, t) and 
written as 


Vst{r,t) = ^ 




r — r' 


Thus the correlation takes the form 


{Vst{r,t)Vst{r',t))s = 


f (Coo(x,t)Coo(x',i)), 


/ 


■dx'dx. 


Using Eqn. Q for the dehnition of ^oo 

{Coo{x)Coo{y))s = Noooo{x,y), 

and Eqn. (fTOl) 

Nooooix,y) = l{{ioo{x),ioo{y)}), 


(17) 


(18) 


(19) 


( 20 ) 


where {.} is the anti-commutator, (.) is the expectation value and too = ^oo ~ {Too), we arrive at 

^2 ^(xk|{Tbo-(Tbo),Tbo-(Tbo)}|'k) 


{Vst{r,t)Vst{r',t))s = ^ j 


-dx'dx. 


X — r X' — r' 


( 21 ) 
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Using the form Too(x) = c^g = mc^|x)(x| for a single particle, we get 


{Vstir,t)Vst{r',t))s = 


2^2 r 


G^m 






lx — r X — r' 


-dx — 2 


( 22 ) 

J |x-r||x'-r'| J 

^ n _ 


We can work out the master equation for the density matrix, following the methods of [401]. The 
evolution of the state vector in (|16p can be written as follows 
d|T(t)) 


dt 




iVr J dk(24'f(k))24(k) + J dk|(k,t)i(k) 


I^W), 


(23) 


where ^(k) = ^ y/l = (G/^vr^h), (^^(k)) = ('I'|24^(k)|'I') and where 


l(k,t) = 


dre 


—ik.r 


dr 


,?oo(r'T) 


r — r' 


(24) 


is the non-white stochastic noise. Using the techniques of |40||, the non-Markovian master equation 
for the evolution of the density matrix corresponding to Eqn. ()23p can be shown to be 


dp 

dt 


[H, p] + iy/j j dk (^{A^)Ap - pA^{A)^ 


7 J dt'dkd'k'D{'k,'k',t,t')^A{k)A{'k',t'— t)p +pA^lk'^t'— t)A^{k)^, 


(25) 


where p = (|'I')('I'|)s. A{k',t' — t) is the interaction picture operator evolved up to time {t' — t). 
The proof for the derivation of this master equation is given in Appendix I. 

It is challenging to demonstrate decoherence from the above master equation, but an alternative 
method applied by Karolyhazy, which we call the phase variance method, comes to our rescue. The 
equivalence of the phase variance method and the Markovian master equation of Diosi’s model has 
been demonstrated by us in Q and it is expected that this equivalence holds in general for 
non-Markovian master equations as well. 

IV. PHASE VARIANCE METHOD AND PROOF OF DECOHERENCE 


In this section we adopt the scheme used by Karolyhazy in his work, to find the decoherence 
effect due to the stochastic potential in the stochastic SN equation (fTbP . In our analysis, we have 
found that 


where 


= -|^V2T(r, t) + U(r, t)T(r, t) 
at 2m 


2 f Gm f ?oo(r',t)^^, 


V (r, t) = —Gm?‘ J 


r — r' 


/ 


-dr' 


r — r' 


(26) 


(27) 
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We now make an important approximation, in order to be able to make progress. In the 
potential V, the exact state 'k will be replaced by the free wave function, denoted as ip, this being 
the solution of the Schrodinger equation with gravitational back-reaction ignored (i.e. 1^ = 0). Thus 
the gravitational effect is being calculated iteratively. In other words, in our further calculations 
we will use 


V{r,t) = —Gm^ J ■ 




r — r' 


■dr' — 


am I 

J |r —r'l 


(28) 


This has the important consequence that the stochastic Schrodinger-Newton equation becomes 
linear, as is the case for the models of Karolyhazy and Diosi. Thus the deterministic nonlinear 
SN equation has now been replaced by a stochastic linear SN equation. The corresponding master 
equation for the density matrix also becomes linear in the density matrix, in this approximation. To 
our understanding, having a linear master equation accompanying a stochastic linear SN equation 
alleviates the problem of superluminal signalling faced by the original deterministic nonlinear SN 
equation. 

Following the argument of Karolyhazy, we can demonstrate that the solution of this linear 
stochastic equation can be given in the form 


T(r, t) = V’(r, t) exp{i(j)st{r, t)) 


(29) 


where '0(r,i) is the free solution without any gravitational back-reaction and is the phase 

of the actual solution with respect to the free wavefunction, and is given by. 


(l)st{r,t) =-^ J V{r,t')dt' 


(30) 


For completeness, we outline the justification of this result of Karolyhazy. The Schrodinger 
equation in the presence of the stochastic potential is given by 




(31) 


where the time dependent Hamiltonian is given by, 


H{r,t) = Ho{r) + V{r,t) 


(32) 


where the second term denotes a small perturbation around the time-independent Hamiltonian. To 
solve this equation, Karolyhazy considers the perturbation to be switched on adiabatically (see, for 
instance, Karolyhazy’s detailed paper of 1990, pg. 223 B). Hence we can aesn.e H.e ..aiabaCic 
approximation” to be valid i.e., if a system starts with a certain state, it will remain in that same 
state after the perturbation is introduced adiabatically. 






Now, let the initial state be 'I'(r, 0) = ^p{r,0) where 'ip(r,t) is the solution without the pertur¬ 
bation so that = 'tp{r,0)exp{—iEt/h). If now, the perturbation is turned on adiabatically, 

the system essentially remains in the same state, with only some phase factors introduced, similar 
to the unperturbed case. Using the adiabatic approximation, we can write (see Griffiths [dll] Chap. 
10, Eqns. 10.12, 10.13 and 10.23) 


'I'(r, t) = 0) exp[i0(r, t)] exp[i 7 (r, t)] 


(33) 


where 6{r,t) = E^t -|- fg V(r,t') dt' , Eq being the unperturbed eigenenergy and V is due to 
perturbation. The phase 7 (r,t) is given by Eqn 10.22 of Griffiths Hence the above equation 
can be rewritten as 


i 1 

T(r, t) ='i/’(r, t) exp —- / V{r,t')dt' exp[i 7 (r,t)] 

L h Jo 




The last phase term 7 (r, t) (called the geometric phase) can be neglected as it depends weakly on 
r (see Karolyhazy jd]). Thus we recover the desired solution written above in Eqns. (1291) and (1301) . 

A few remarks are in order, with regard to the application of this phase variance method to 
the present case of the stochastic SN equation (IMj) . Both the terms in the potential (|28]) . namely 


the mean as well as the fluctuation, are treated on the same footing 


32l ]. as in principle there 


can be circumstances when they are comparable. Thus, from the viewpoint of application of the 
adiabatic approximation assumed while writing down the ansatz (|29p it is assumed that the total 
potential can be treated as a perturbation on the free particle. While this is definitely evident for 
microscopic objects, the assumption maybe considered to be a reasonable one in the macroscopic 
case as well, so long as the gravitational field in question is sufficiently weak. The plausibility of 
the results that we derive below serves to justify the validity of the ansatz ([29]). 

The significance of the form (I29p of the solution to the stochastic SN equation is the following. 
The phase at any spatial point r is a stochastic variable. If we consider two fixed spatial 

points ri and r 2 , the phase difference between the quantum states at these two points is also 
stochastic, and the variance of this difference grows with time. When the variance becomes of the 
order vr^, we say that the states at these two points decohere. This gives the decoherence time 
as a function of the spatial distance |ri — r 2 |. This decoherence criterion has been shown to be 
equivalent [l^ to the one defined using a Markovian master equation, as in Diosi’s model, and we 
assume that the equivalence holds in the present non-Markovian case as well. The plausibility of 
the assumption is supported by the nature of the results we find below, which tally with results 
obtained by others earlier. 
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We now calculate the phase variance between two points ri and r 2 with the formula 


^ - (j)stir2,t)f)s 


(35) 


We find the time t for which the above quantity would be ~ vr^ and that will give us the decoherence 
time. Calculating the variance using the definitions (1301) and (f28]l we get, 


2^4: ^ 


. ,2 


r',t')nV’(r",niV3 


|ri — r'l |ri — r"| 


d^r'd^r''dt'dt” 


+ 


II 




|r2 - r'||r2 - r" 




|ri — r'l |r2 — r'' 




+ 


G^m 


V 


(r•^tOP 

,/|2 


- 2 


/ 


|ri — r 




d^r'dt'dt" + 


■d^r'dt'dt" 


IV^(r•^tOP 

|r 2 — r'P 


■d^r'dt'dt” 


ri - r'||r2 - r' 


(36) 


where the time integration is done up to a time T and the volume integrations are over all of space. 

To calculate the phase variance, we use the free gaussian solution of the wave function for a 
point particle of mass m 2^ which is given by [we recall from our remark above (|28p that we are 


finding the gravitational effect perturbatively] 


t/;(r,t) = (vra^) ^ g^p [ 

\ ma^ / \ 


2 a 2 (l + 

^ ma^ ' 

Substitution of this free wave function in the above variance formula (|36p gives the result 


(37) 


A(/.^ = 


2 _ rj,2^ 1 


^2 


_ J _\2 

n r2' 


(38) 


Details of this derivation are given in Appendix II. We demand the phase variance to be of the 
order 7r2 1. This gives. 


h 


1 


Gm? 


I j_Ll 

In r 2 I 


(39) 


Hence 


T = 


h 


AE = 


Gm? Gm? 


(40) 


|A£'|’ ri r 2 

Thus, from the stochastic Schrodinger-Newton equation we have deduced a result similar to the 
Diosi-Penrose criterion, namely that the decoherence time is of the order of the inverse of the grav¬ 
itational potential energy difference at the two locations. This strongly suggests that gravitational 
decoherence drives the system to one of the solutions of the SN equation, as originally proposed 
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by Diosi and Penrose. It is interesting that this result has been derived starting from the Einstein- 
Langevin equation, and it highlights the significance of considering quantum fluctuations in the 
stress tensor, while also encouraging us to believe that stochastic gravity is intimately connected 
with the stochastic SN equation. 

Here, it is evident that we certainly cannot talk of the coherence length as simply the distance 

— ^ 2 ! between the two points because the gravitational potential energy difference is involved. So, 
unlike the other gravity induced decoherence models, where the decoherence time and coherence 
length depend upon the distance between two points, here they actually depend upon the difference 
1^^-3-|. So for this case, we introduce a length scale Oc, which is the coherence length, defined as: 


1 


Then from the above calculation we can write. 


1 

n 


1 

r2 


T = 


httc 

Girfi 


Again we know that T ~ ma^/ti. These two expressions together give, 


dn — 


Gm? 


(41) 


(42) 


which is a well-known result. 

In our calculation, we have found that the decoherence time between states at the two points 
ri and r 2 depends on the difference between the gravitational potential energy at these two points. 
This is somewhat similar, although not exactly same as, the so called Diosi - Penrose criterion. 
The latter says that if we have two stationary states represented by two mass configurations then 
their superposition will decay to one of the solutions of the SN equation in a time scale r Ri h/ A.Eg 
where AEg is the gravitational self-energy of the difference between the mass distributions of the 


two states 


0,£i. 


Penrose suggests in 


ITi ] that if we have two mass distributions representing 


A = -47rG 


the two states characterised by the mass densities g and g' respectively, then the measure of 
incompatibility between them would be 

(^.(x) - ^^'(x))(^(y) - g'{y))/\^ - y| £y (43) 

He suggests that the decay time would be of the form r ~ h/ AEg where AEg is A or some multiple 
of that quantity. Diosi also obtains the same form, as has been given in [14l . |42 |. where he finds 
that the decay time should be r ~ h/AEg, with AEg having a form 

[/(x|^) -/(x|A')][/(x'|A) -/(x'|A')]/|x-x'|d^xd^x' (44) 

Here /(x|A) denotes the mass density at x for a configuration denoted by X. 


AEg = G 
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V. SOME QUALITATIVE ESTIMATES FROM THE SN EQUATION 


In an effort to obtain some estimates as to how the coherence length Oc depends on the mass m 
and size R in case of an extended object, and in order to compare with the results of Karolyhazy, 
and of Diosi, we make some qualitative estimates. 

Let us recall first the case of the point particle treated in 


22l | according to the SN equation: 






r — r' 


dr')'i/)(r, t). 


(45) 


An estimate of the critical length and decay time can be obtained from this equation considering 
the self-gravitation of an expanding wave packet. If we start with a spherically symmetric Gaussian 
wave packet of width a, 

0 ) = 


exp 


r 


(46) 


then after time t, the wave function evolves, via Schrodinger equation (in the absence of any 
gravitational potential), as the state given by (1371) . The radial probability density is at a maximum 
at rp = a^J 1 -|- So, the peak shifts with time with an acceleration given by the following 

equation: rp = h?/m?rp. This gives the acceleration of the wave packet when it expands freely. 
Now, the acceleration due to gravity for a point mass m at a distance Vp is rp = Gm/rp. Let, 
after a certain time, these two accelerations become equal. We call the width of the wave packet 
in such an equilibrium as aequU- Then equating the two accelerations we get, for a mass m, 
m = /Gaequii^ ■ So, for a given mass, one can calculate this equilibrium width. If rp < a^quii 

then usual quantum evolution dominates while for > a^quU gravity becomes more significant and 
collapse of the wave function takes place. Hence, aequU is the same as as the critical length ac- 
We now focus our attention on the behaviour of the wave function at t = 0. At t = 0 we have 
rp = a. Now if we start with a mass m such that /m?a^ > then the Scrodinger-like expansion 
keeps on accelerating unless the two effects become equal. In that case we have Uc > o. But if 




< 


Gm 


then the wave packet starts contracting right from the beginning and we get ac < a. 


This can be summarized as below: 

1. If m < ( ^ I then ac > a 

then ac < a 

So, for an already contracted wave function, we must have in > (^) 
threshold mass for collapse. 



,1/3 

Ga^ 

) 

'' 

1/3 

Ga . 

1 


m = 


Ga 


1/3 


is the 
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Critical Length for Extended Objects : 

We have already seen that the critical length for a single point-like mass m is given by Uc = 
h^/Gm^ Now we calculate the same for an extended object of mass m and size R. We consider 
two cases: (i) Oc S> R- In this case, the gravitational acceleration on the wave packet which 
extends outside the mass m is, g = Gm/Vp. For the critical length Oc we have, /m?a^ = Grajc?^. 
This essentially gives the same critical length as in case of a single point mass Oc = f? jGm'. (ii) 
Oc ^ R. When Oc <C R, the wave packet lies inside the extended mass m. Considering the mass 
m having uniform density within the radius R we get the gravitational acceleration on the wave 
packet as, g = Gmrp/R?. Again at the critical length, = GmadR?■ This gives the 

critical length as: Uc = h/Gm?Y^^R^^^■ Interestingly, this same expression for Oc is obtained by 
Diosi in his treatment of gravitational decoherence, while it differs from the expression obtained 
by Karolyhazy. What this is telling about possible similarity between the stochastic SN equation 
and the Diosi approach is not clear to us at present. 

Putting Oc ~ i? gives the transition between micro region and macro region: 

• micro region(m^i? <C |t): critical length Oc = 



• Transition occurs at i? ~ 


At rp = Qc the two opposite accelerations become equal and so, the wave function evolves at a 

/ .2 \ 1/3 

constant rate. Let us consider a mass for which m is below the threshold value i.e, m ) 

Initially when Vp < Uc, the wave packet expands at an accelerating rate until it reaches rp = Uc- At 
this point, as the two forces become equal, it stops accelerating and starts expanding at a constant 
rate. As soon as rp becomes greater than Oc, again, gravity dominates and the wave packet shrinks 
to size rp = ttc- Again at this point onwards, it keeps on contracting at a constant rate until < Oc 


where again Schrodinger like expansion takes over. Thus, there should be an oscil 
the equilibrium width Oc which has been found through numerical simulations in 


M- 


ation around 


VI. CONCLUDING REMARKS 

We believe that the stochastic SN equation we have found makes a useful contribution to the 
subject of gravity induced decoherence. The equation is possibly robust, because its origin lies in 
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the well-defined theory of stochastic semiclassical gravity. It correctly predicts the Diosi-Penrose 
decoherence criterion which relates decoherence time to the gravitational potential energy difference 
between the two points under consideration. The equation is on the same footing as the stochastic 
Schrodinger equations of Karolyhazy and Diosi, especially in its linearised form obtained after 
using the approximate stochastic potential (I28p . which choice makes the equation linear. The key 
difference of our equation is that it deals with the stochastic effects of self-gravity, whereas the 
Karolyhazy / Diosi models deal with stochastic effects induced by intrinsic space-time uncertainty. 
How these two effects relate to each other, and whether a more general treatment should include 
both effects remains a subject of future investigation. Like the master equations in the CSL and 
Diosi models, which resemble collisional decoherence models j3 |. our master equation in principle 
does the same, although explicit details remain to be worked out. 

Of course the full correlation function (l2Tp in our case is state dependent, unlike in the case of 
the stochastic potentials introduced by Karolyhazy and by Diosi. Furthermore, the noise is not 
white, unlike in the case of Diosi’s model, where the noise is white in time. And yet, the final 
physical results match in the two cases - the reasons for this intriguing feature are not clear to us 
at present. 


In an interesting recent paper Nimmrichter and Hornberger 


3.‘ll | have also proposed stochastic 


extensions of the (regularised) SN equation. One of their goals is to make the SN equation com¬ 
patible with the no-signalling requirement. They suggest two possible stochastic extensions: one 
via discrete jumps described as a Poisson random process, and another as a continuous white noise 
(Wiener process). It is shown that nonlinearity, and hence superluminality, is avoided at the level 
of the stochastic average described by the density operator, which obeys a linear master equation. 
Our work differs from theirs in the key aspect that we derive the stochastic noise starting from the 
Einstein-Langevin equation, and this noise is different from the ones they propose. However, as 
noted above, if we work with the approximate potential (I28p [gravitational back-reaction is calcu¬ 
lated iteratively] we also obtain a linear master equation, consistent with no signalling. Also, the 
noise that we introduce seems to lead naturally to the gravitational decoherence effect described 
by Diosi and Penrose. 

The A^-particle stochastic SN equation can be constructed following the method described in 
(their Eqns. 6-12) and by using for the potential the form ()15p given above. Generalizing the 
decoherence result obtained above, we expect this equation to produce localisation of macroscopic 
masses. Thus the two-particle stochastic SN equation can be expected to localise the two particles 
and correctly reproduce their classical Newtonian gravitational interaction. 
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An aspect which we have not yet investigated is the magnitude of stochastic heating brought 
about by the stochastic terms present in the new equation. It is known in the Diosi model that 
there are divergences in the master equation in the point particle limit, resulting in a divergent 
heating rate in this limit. One possible w ay to avoid this divergence is to introduce a cut-off (see for 


instance the detailed recent discussion in 


30 ] ). It is perhaps possible that this divergence problem 


does not arise in our model because the non-white noise in the non-Markovian equation induces a 
finite correlation length; however at the moment this is a conjecture, which needs to be scrutinised. 

Further investigations should include studying the stochastic SN equation for states other than 
Gaussian wave-packets, and for extended objects. It will be useful also to see if there is some 
way to convert the stochastic equation into a collapse equation, and to compare it with CSL. One 
should also attempt to work out how the stochastic SN equation can be tested in experiments such 
as molecular interferometry and via gravitationally induced random diffusion \m. 45|. 
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Appendix I : Master equation for the density matrix 

We have from Eqn. (I23p 

d|vk(t)) 


dt 


+ / dk(it(k))A(k)+iV7 / dk|(k,t)A(k) 




(47) 


It implies that 

d\'i>) = -^Ho\'i>)dt + i^ j dk{A\k))A{k)\'f)dt + i^ j (ikf(k,t)i(k)|T(t))dt, 

and, on taking the conjugate of the equation 


(48) 


d(^| j dk(i(k))i^(k )j dkr(k, t)i^(k). (49) 

The stochastic SN equation preserves the norm of the statevector |'k). To see this we compute 
the differential of the stochastic quantity ('k|'I'): 


d{'k\'k) = d(T||T) -b (^|d|^) 


(50) 
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From equations (1481) and (|49p we have 

j d\i{A)A^\^) - dt{-^\i^ j dkiA^^) 
- m^Ho\'^)dt + i^{^\ j dk{A^)A\^)dt + i^{^\ j dkCA\'k)dt 

= i^ J dk(^|(^|i|^)dt 
= 0 

Let Pst = |^)('h|. Then 

dpst = d|^')(^'| + |^)d(^| 

= -^[Ho,Pst]dt + iy/^ j dk (^{A'^)Apst - PstA^A)^ dt 

+ iy/7 J (^Mpst - PsttA^^ dt. 

The density matrix of the system is the stochastic average of the pst 

P = {Pst)s- 


(51) 


(52) 


(53) 


We now take the stochastic average of Eqn. (|52p and retain terms upto the order of ^/y. It simply 
means that (A) Ri (TolAlTo) where we have the expansion IT) = ITq) + v^lV^i) + 7 l'I' 2 ) + • • • 


{dpst)s = \ - -^[Hq, pst]dt + iyjPy / dk 


st Pst 


dt 


+ *V7 j dk {A^pst - PstCA^^ ( 

^ =-^[id,p]+i^/7 / dk(^{A^)Ap-pA^A)'j+iy/j I dk(^A(^ipst'j(^Pstt^j^A^'j 


According to the Furutsu-Novikov formula 


F[|]f(k,t)^ = j dk' J dt'D{k,k',t,t') 


5|(k',d)/ ’ 


(54) 


(55) 


where F[^] is an arbitrary functional of the gaussian stochastic noise and D(k,k' ,t, s) is the cor¬ 
relation function of ^(k, t) 


D{k,k',t,t') = {^{Kt)A{k',t'))c 

= (^kk' J dridr2e-*^'^ie-*^''^2 J ,t) Coo{r",t') 

= kk' j dridr2e-*^'^ie-*^'-’'2 j 


ri — r'l |r2 — r"| 
./ „// + +/\ 


( 56 ) 




|ri — r'l |r2 — r"| 
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where Afoooo(^^ r", i') is given by Eqn. Q. 

Let F[^] = psf We assume a system with spherical symmetry. Then, 

= A-k J dk'k''^ J dt'D{k,k',t,t') . (57) 

We now transit to the interaction picture where the operators and state vectors are defined as 
i(k,t) = \A>{t))i = F^\A>{t)). (58) 


We, then, expand IT)/ with the parameter as follows 

|T)/ = |To)/ + + 71 ^ 2 )/ + • • • (59) 

In the interaction picture Eqn. (1471) takes the form 

= iVi/rfk (atk)) + «k, *)) A(k, ()I'I'(*))/- ( 60 ) 

On plugging Eqn. (I59p in Eqn. (IBOp and collecting the terms of order 0, and 7 we get 

0 ( 0 ): 


— I 

dt' 


To)/ = 0 


(61) 


0 (^ 7 ): 


^|Ti)/ = z J dk (^(i^(k))+|(k,t)^ i(k,t)|To)/ 
l^i(*))=* J dt' j dk (^(i'f(k))+^(k,s)^ i(k,t'- t)|To(t)) 


(62) 


0 ( 7 ) : and so on... 

The density matrix in Schrodinger picture pst = |'h)(T|, when expanded perturbatively upto 
order is 


Pst = |To)(To| + V7(l^i)(^o| + l^o)(^i|)- 


(63) 


Using Eqn. (| 6 ^ . we can evaluate 

and 


dpst 

df*(k', t') 


,t' - t)), 


(64) 


(65) 
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where = |'I'o)(^'o|- 

Hence from Eqn. (153), we have 


= Airiy/j j dk'k''^ J dt'D{k,k',t' — t)p] 
= Airi-s/j J dk'k''^ J dt'D{k,k',t,t'){A(k',t' — t)p) 


( 66 ) 


Finally from Eqn. (|54p we obtain 


^ +iV7y dk(^{A^)Ap-pA^A)'^ J dk(^A(^ipst'^^-(^Psti'^^A^'^ 
-■-^[H,p]+i^/A J dk(^{A^)Ap- pA^A)^ 

'Y J j dkdk'D{k,k',t,t')^A{k)A{k',t'— t)p + pA^{k',t'— t)A^{k) 


(67) 


where D{k, k',t,t') = lQ'K‘^k‘^k'‘^D{k,k 'We note that the calculation has been carried out 
perturbatively, by expanding the state in powers of and is exact up to order y/ 7 . 

Appendix II: Calculation of phase variance for gaussian state 

Let us consider the first term in ()36l) . The integration can be separated into two integrals of r' 
and r”. The integration can be performed in spherical polar co-ordinates. Similarly, all the other 
terms can also be calculated. Some steps of the calculations are as follows. The hrst term can be 
separated into two integrals of r' and r" of the form f f d^r'dt'. We have, 




Now let ri — r' = r so that r'^ = rf + — 2rir cos 9. This gives 


,.'2 


rriFa^ ^ 


r', t')\^ = TT exp { — ^{r‘i + r'^ — 2 rir cos 0 )) 


( 68 ) 


(69) 


a-^ \ 

where a = 1/(1 + ^). 

The integration is now straightforward. Other terms can also be calculated in the same fashion. 
Small time approximation: In our calculation, in order to obtain an analytical result, the time 
integration has been calculated using small time approximation i.e. the interval T is such that 


_„2/_2 h^ri r/ , N 

5 1 / ^ < erf n/a) 


(70) 


where the error function erf(x) is given, along with the imaginary error function erh(x), used below, 
as 


erf(a;) = 


A r 

Jo 


e ^ dt 


erh(x) = —ieii{ix) = —j= / e* dt 

'TT Jo 


(71) 
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Th condition given in (17011 is well satisfied for ri/a —)• oo. 
After using small time approximation we get 


A(/.2 = 


+ 


3G2 


m 


G'^m^ 


rj-iZ 


2 , fert (n'lV + 1 fert (Ilif - in') erf (U 


4R2 


TT 


a 7 7 


n 


a / J 


nr 2 


a 7 


\ a 


ria 


1 -e 7^ erfi ( — ) H-e 7^ erfi ( — 


TT 


r2a 


r2 


( 72 ) 




The last term in the above equation has not been calculated explicitly as it was very complicated. 
So we used the argument that since it was initially obtained by breaking a squared term, the final 
result should be a perfect square. Putting the last term in this form gives the desired result. 

We now calculate term by term in Eqn. (1721) in the limit ri/a —)• oo and r 2 /a —>■ oo. First we 
look at the terms involving ri. So, the hrst and fourth terms of the above expression together give, 

r 2 - 


4^2 


r2 


0f!le-'-?/“'erfif—) +3eii^(— 


As — ^ oo, the term inside 


goes to 4. So, we have, this term as 

r 2 


4^2 r 


•4 


Similarly, by adding second and fifth terms, we get, 

G^m^ r2 
4^2 TI" 


•4 


(73) 


(74) 


(75) 


Now let us consider the third and sixth terms together: 


G^m^ 2r2 r 


3erf 


(V)-(V) - (.0) 


Ah? rir2 L 

Again the term in square brackets goes to 4 for ^ oo and ^ n- oo. We get, 

•4 


G2m4 2 r 2 


4h2 


rir2 


(77) 


Finally, adding all the terms gives the desired result 


2 G‘^rn^ m2 


ri r2 7 


(78) 


which is Fqn. (f38l) . 

The calculations can be repeated for the case ^ ~ 1 and ^ ~ 1 and the results come out to 
be the same. The only problem in this limit seems to be that the small time approximation is not 
valid for microscopic masses. For large masses, it is still valid. 


19 





















REFERENCES 


[1] F. Karolyhazy, Nuovo Cimento A 42, 390 (1966). 

[2] F. Karolyhazy, A. Frenkel, and B. Lukacs, in Quantum concepts in space and time, edited by R. Penrose 
and C. J. Isham (Clarendon, Oxford, 1986). 

[3] F. Karolyhazy, Magyar Fizikai Folyoirat 22, 23 (1974). 

[4] F. Karolyhazy, in Sixty-two years of uncertainty, edited by A. Miller (Plenum, New York, 1990). 

[5] F. Karolyhazy, in Fundamental problems of quantum physics, edited by M. Ferrero and A. van der 
Merwe (Kluwer Acad. Publ., Netherlands, 1995). 

[6] F. Karolyhazy, A. Frenkel, and B. Lukacs, in Physics as natural philosophy, edited by A. Shimony and 
H. Feshbach (MIT Press, Cambridge, 1982). 

[7] A. Frenkel, in Quantum mechanics - a half century later, edited by J. L. Lopes and M. Paty (Reidel, 
Dordrecht, 1977). 

[8] A. Frenkel, Found. Phys. 20, 159 (1990). 

[9] A. Frenkel, in Fundamental problems of quantum physics, edited by M. Ferrero and A. van der Merwe 
(Kluwer Acad. Publ., Netherlands, 1995). 

[10] A. Frenkel, Found. Phys. 32, 751 (2002). 

[11] A. Frenkel, in Experimental Metaphysics, Vol. 1, edited by A. Shimony, R. S. Cohen, M. Horne, and 
J. Stachel (Kluwer Acad. Publ., Boston, 1997) p. 39. 

[12] L. Diosi, Physics Letters A 120, 377 (1987) 

[13] L. Diosi, J. Phys. A 40, 2989 (2007). 

[14] L. Diosi, Phys. Rev. A 40, 1165 (1989) 

[15] L. Diosi and B. Lukacs, Annalen der Physik 499, 488 (1987). 

[16] S. Bera, S. Donadi, K. Lochan, and T. P. Singh, Preprint , aXiv: 1408.1194 [Found. Phys. (to be 
published)]. 

[17] R. Penrose, Gen. Rel. Grav. 28, 581 (1996). 

[18] R. Penrose, Phil. Trans. R. Soc. London A 356, 1927 (1998). 

[19] R. Penrose, in Mathematical physics 2000, edited by A. Fokas, A. Grigoryan, T. Kibble, and B. Ze- 
garlinski (Imperial Gollege Press, London, 2000) p. 266. 

[20] L. Diosi, Phys. Lett. A 105, 199 (1984). 

[21] D. H. Bernstein, E. Giladi, and K. W. R. Jones, Mod. Phy. Letts. A 13, 2327 (1998). 

[22] D. Giulini and A. Grofiardt, Glass. Quant. Grav. 28, 195026 (2011). 

[23] R. Harrison, 1. Moroz, and K. P. Tod, Nonlinearity 16, 101 (2003). 

[24] R. Harrison, 1. Moroz, and K. P. Tod, Class. Quant. Grav. 15, 2733 (1998). 

[25] R. Ruffini and S. Bonazzola, Phys. Rev. 187, 1767 (1969). 


20 





[26] D. Giulini and A. Grossardt, Glass. Quant. Grav. 29, 215010 (2012). 

[27] D. Giulini and A. Grossardt, Glass. Quant. Grav. 30, 155018 (2013). 

[28] B.-L. Hu, J. Phys. Gonf. Ser. 504, 012021 (2014). 

[29] G. Anastapoulos and B. L. Hu, New J. Phys. 16, 085007 (2014). 

[30] M. Bahrami, A. Grossardt, S. Donadi, and A. Bassi, New J. Phys. 16, 115007 (2014). 

[31] S. Colin, T. Durt, and R. Willox, Class. Quant. Grav. 31, 245003 (2014). 

[32] B. L. Hu and E. Verdaguer, Living Rev. Relativity 7, 3 (2004). 

[33] S. Nimmrichter and K. Hornberger, Phys. Rev. D 91, 024016 (2015). 

[34] B. L. Hu and C. Anastapoulos, Preprint , arXiv:1504.03103 (2015 Probing a gravitational cat state). 

[35] A. Bassi, K. Lochan, S. Satin, T. P. Singh, and H. Ulbricht, Rev. Mod. Phys. 85, 471 (2013). 

[36] S. Gao, Stud. Hist. Philos. Mod. Phys. 44, 148 (2013). 

[37] T. P. Singh, J. Phys. Gonf. Ser. 626, 012009 (2015). 

[38] B. L. Hu and E. Verdaguer, Class.Quant.Grav. 20, R1 (2003). 

[39] R. Martin and E. Verdaguer, Phys.Lett.B 465, 113 (1999). 

[40] S. L. Adler and A. Bassi, Journal of Physics. A. Mathematical and Theoretical 40, 15083 (2007). 

[41] D. J. Griffiths, Introduction to quantum mechanics, 2nd ed. (Pearson Prentice Hall, New Jersey, 2014). 

[42] L. Diosi, Braz. J. Phys. 35, 260 (2005). 

[43] B. Vacchini, J. Phys. A 40, 2463 (2007). 

[44] B. Collett and P. Pearle, Found. Phys. 33, 1495 (2003). 

[45] S. Bera, B. Motwani, T. P. Singh, and H. Ulbricht, Scientific Reports 5, 7664 (2015). 


21 


